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A solid - amorphous or crystalline - is defined by a finite shear modulus while a fluid lacks such. 

We thus experimentally investigate the elastic properties of a colloidal glass former near the glass 
transition: spectroscopy of vibrational excitations yields the dispersion relations of longitudinal and 
transverse phonons in the glassy state. From the long wavelength limit of the dispersion relation 
we extract the bulk and the shear modulus. As expected, the latter disappear in a fluid and we 
measure a clearly resolved discontinuous behaviour of the elastic moduli at the glass transition. 

This not only determines the transition temperature To of the system but also directly addresses 
recent discussions about elasticity during vitrification. We show that low frequency excitations 
in our system are plane waves such that continuum elasticity theory can be used to describe the 
macroscopic behaviour. 


In general, there are several ways to define glassy sys¬ 
tems mm- For instance, a characteristic change in the 
thermodynamic properties volume and enthalpy upon 
cooling a fluid may be used to define the glass transition 
temperature Tq. Another property often investigated 
for molecular or atomic glasses is the viscosity t]. Taking 
into account the rapid slowing down of dynamics, the 
glassy state is reached when the viscosity exceeds 10^^ 
Pa-s upon crossing the freezing temperature. However, 
all of these definitions are difficult to apply to soft matter 
systems. In studies on the mesoscopic or microscopic 
scale, access to the macroscopic property r] often is 
difficult, if not impossible. Secondly, these systems are 
“softer” up to a factor of 10“^^ (2D) or 10“^® (3D) which 
rules out the idea of an universal, viscosity-dependent 
categorization. Thus, other ways have been suggested 
to pin down Tq [ 1 ]. A manifestation of solidity is the 
emergence of rigidity, expressed by the shear modulus 
/r becoming finite. Microscopically, prerequisites for 
solidity are long ranged dynamical correlations. While it 
is easy to understand how such long range correlations 
emerge in a crystal (based on the concept of long ranged 
periodic order), solidity is much less easy to grasp on 
an amorphous, disordered background. As the shear 
modulus is zero in fluids, we expect a significant change 
of /r at the onset of vitrification mm- Although ^ is 
a macroscopic property, it is measurable locally via the 
equipartition theorem which makes this concept well 
suited for soft matter systems mm- From an experimen¬ 
tal point of view, this approach is advantageous because 
it relies on measuring the change of a quantity from zero 
to a finite value. On the other hand, the usual approach 
of measuring the divergence of dynamical quantities like 
the viscosity is much harder to realize experimentally be¬ 
cause it requires access to a much larger parameter space. 

Elasticity of complex fluids and their viscoelastic 
behaviour has been tackled mostly by rheological mea¬ 
surements so far [7]. In general, these methods allow 
the study of frequency dependent storage G'{uj) and 


loss G”{uj) modulus. The former describes the elastic 
contribution, while the latter covers dissipative processes 
upon shear. Thus, a non-vanishing G'{uj) > G”{uj) in the 
limit w —>■ 0 indicates a solid. These quantities have been 
investigated extensively by experiment as well as theory 
El and from the data, a clear distinction between 
solid and fluid state might be drawn. This is also 
true for flow curves, which might be not only acquired 
by conventional rheometers but also by microfluidic 
rheology m- Here we address solely the elasticity by 
internal fluctuations without external driving. 

While it has been shown recently that the zero- 
frequency shear modulus lim,^_,,o G"(a;) = /i is a good 
quantity to observe the glass transition m, there has 
been an ongoing debate on the behaviour of /r when an 
amorphous solid forms. Replica theory approaches [12] 
predict a continuous growth from zero upon cooling. 
This agrees with the view from granular systems, where 
theory predicts critical fluctuations close to the jamming 
transition of these athermal systems which also lead to 
an algebraic growth of /x from zero m- On the other 
hand, another recent replica theory calculation reveals 
a discontinuous jump at the glass transition [T4lH6j . 
This complies with mode coupling theory predictions 
m where ^ is connected to a non-ergodicity parameter 
appearing in a stress auto-correlation function. A 
discussion about the discontinuous elasticity in the 
context of Random First-Order Transition is given in 
[IB] . In addition, simulations of hard discs also showed a 
discontinuous onset of the shear modulus El by access¬ 
ing the system with the equipartition theorem. Similar 
results from novel four-point correlation calculations 
back this up further [19]. 

Here, we evaluate microscopy data from a two dimen¬ 
sional colloidal model system with soft interactions in 
equilibrium. This gives us access to information on the 
single particle scale as opposed to atomic or molecular 
systems, where only ensemble data may be studied. By 
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determining the displacement field and subsequent anal¬ 
ysis in Fourier space, the dispersion relations of acoustic¬ 
like excitations are obtained and provide insight into the 
macroscopic elastic properties of the system. As the sys¬ 
tem is cooled, we are able to confirm the discontinuous 
nature of the shear modulus at the glass transition. 


METHODS 


Video microscopy provides access to the full phase 
space information. We record the trajectories of about 
2300 colloidal particles confined to two dimensions at a 
flat water-air interface. The species A (diameter a = 
4.5/im) and B (cr = 2.8 fxm) have a relative concentra¬ 
tion of ^ = Nb/{Na + Nb) « 50% where and Nb 
are the number of particles of both species in the field of 
view. This prevents the system from crystallization. 

The superparamagnetic nature of the particles lets us 
control the particle interactions in situ by an external 
magnetic field H. This is expressed by the dimensionless 
system parameter 


/To 

47r kBT 


{^■XB + {l-0-XAf , 


( 1 ) 


which effectively acts as an inverse temperature. Here, n 
denotes the area density and is computed via a Voronoi 
tessellation. xa,b represent the susceptibilities of species 
A and B, respectively. After equilibration at low T, the 
system was cooled down stepwise. For different states up 
to 3.8 • 10® snapshots were analyzed. With a frame rate 
of « 2s“^, sampling times up to some 10® seconds were 
achieved. This is sufficiently long to probe dynamics even 
in the relaxation time regime Tq, for highest system pa¬ 
rameters F. Additional details of the setup are described 
elsewhere [301 [H]. 

We follow the well established data evaluation scheme 
described in [22]. The harmonic part of the systems po¬ 
tential energy is bilinear in the displacements u times the 
“spring constant” given by [32 


U . ( 2 ) 

9,a,/3 


For the monolayer the dynamical matrix (q) is a 2 x 2 
matrix and describes the elastic coupling between parti¬ 
cles. It is given by a sum of the second derivatives of the 
interparticle potential and a, /3 S 1,2 index the two spa¬ 
tial dimensions. Elasticity naturally originates from the 
spring constants of the system. These are the eigenvalues 
Ss of the dynamical matrix we are interested in. 

The displacement Ui(t) = — of particle i is taken 

relative to its equilibrium position and its Fourier compo¬ 
nent is given by u{q) = We calculate 

those positions to be the center of mass of the trajectory 
{fi) = 1/M • nit), where the average is taken over a 


time interval At significantly longer than any short time 
relaxation process. In the equation, M corresponds to 
the number of snapshots during the time interval At and 
in an amorphous solid those time averaged positions con¬ 
verge to the quasi-equilibrium positions. The expression 
for (f/) holds as long as At is smaller than the relaxation 
time Ta and particles are confined in their “cage”. was 
extracted from MSD measurements for high F. 

In an amorphous system we cannot calculate the dynam¬ 
ical matrix a priori, but using the equipartition 

theorem we can measure it’s eigenvalues <5^- The equipar¬ 
tition theorem reads {U) = N ■ kBT where N is the num¬ 
ber of particles, each with two translational degrees of 
freedom. With Eqn. it follows 

kBT = 6s{^ ■ (3) 

The subscript s denotes polarizations ||,T and depends 
on the relative orientation of u and q. For q \\ u and 
q T u we study longitudinal and transverse behaviour, 
respectively. 5a{(^ is termed “dispersion relation” 
but we keep in mind that in our system, oscillatory 
motion is overdamped due to the surrounding solvent. 
By applying this method, we implicitly assume the 
investigated modes to be plane waves. This holds in 
the long wavelength limit corresponding to classical 
continuum elasticity theory [33] ■ There, the 2D system 
is considered to be a homogeneous structureless solid. 
Theoretical background on this approach can be found 
in[TI|. 


From the dispersion relations, the elastic constants can 
be derived in the long wavelength limit [Sj. In a two 
dimensional isotropic solid, the elasticity tensor C^^aT 
possesses only two independent elements. These two el¬ 
ements can be expressed via the Lame coefficients /r and 
A of continuum elasticity theory. We extract the Lame 
coefficients from the “dispersion relation” as: 


kBT 

Oq (2/i -|- A) 


kBT 

alfj. 


liin[4(|M||(^)n] , 

q-^o q‘^ 

lini[4(|Mu(^)P)] • 

q‘‘ 


(4) 

(5) 


Here, ag is the mean interparticle distance, defining the 
area density as n = I/oq. In two dimensions, ^ gives 
the shear modulus and fi + X = B the bulk modulus. 
As we extract the elastic constants for (f —>■ 0, we are 
only accounting for long wavelength modes in elastic 
continuum theory, which justihes our approach even in 
disordered materials. 


Recent work on vibrational properties of colloidal 
glasses relied on the diagonalisation of the covariance 
matrix of the displacement field [351 - [37| . From the spec¬ 
trum of eigenvalues, this two dimensional principal com¬ 
ponent approach derives the density of vibrational states. 
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FIG. 1. Mean square displacement (msd) for both species of 
particles at different temperatures. The kink at high temper¬ 
atures slowly develops into a plateau. However, it is impossi¬ 
ble to distinguish unambiguously between solid and fluid state 
from this dynamical data. The black solid line corresponds to 
free diffusion and the vertical dashed line to the sampling time 
At used in the following analysis to determine the elasticity. 

In all cases (hard spheres and soft PNIPAM particles) 
an anomalous low frequency behavior was found, which 
was connected to the boson peak known from scattering 
experiments in molecular or polymeric glasses |281 l2Hj . 
This excess of modes was attributed to the non-isotropy 
of disordered systems at short to intermediate length- 
scales. In atomic and molecular glasses, recent progress 
showed that long wavelength modes undergo a crossover 
from propagation to localization due to diffraction upon 
decreasing the wavelength [30H34j. By analyzing the 
Eigenvector structure of the covariance matrix, the spa¬ 
tial features of the corresponding modes could be visual¬ 
ized in the colloidal case P5H?7] . It was concluded that 
the low frequency spectrum consisted mainly of local¬ 
ized “swirling” type of modes. Since this does not go 
along with our approach of plane waves in a continuous 
medium, we study the structure of modes via principal 
component analysis as well at the end of the manuscript. 
From the particle displacements u(t), we compute the 
covariance matrix 

Cij = {u,{t)uj{t))t, (6) 

where i and j correspond to all components of all par¬ 
ticle displacements. By diagonalisation we obtain the 
spectrum of eigenvalues c*. Following [2M?7] we use the 
inverse as frequencies of the modes in our systems, pro¬ 
viding a mapping = I/c*. The corresponding Eigen¬ 
vectors reveal the structure of those modes we are inter¬ 
ested in. In the calculations, we follow the procedures 
described in [35]. Because principal component analy¬ 
sis requires large statistics, we use 1.24 • 10® frames for 
constructing the covariance matrix. At this number the 
results have statistically converged. 


RESULTS 

To start off, we present the measured 
mean squared displacement (msd) (Ar^(T)) = 

-I-r) — of the particles in Fig. 

for different system temperatures. Here, N presents 
the total number of particles which are located at 
positions fi at time t. It clearly shows the expected 
slow down of dynamics. At high temperatures, the data 
show a mere kink which slowly develops into a plateau 
for intermediate lag times r. For longer times r^, the 
dynamics exhibits (sub) diffusive behavior corresponding 
to alpha relaxation processes. This becomes clear by 
comparing to the solid black line with slope one, marking 
free diffusion. It is noteworthy that from this dynamic 
data alone, the exact temperature of a transition from 
the fluid state to a solid (amorphous) state can hardly 
be identified. 

Therefore we focus on spectroscopy of acoustic excita¬ 
tions. Fig. [^shows the measured normalized “dispersion 
relations” for different interaction parameters analyzed 
from trajectories of big particles. Filled and empty 
symbols represent spring constants for longitudinal and 
transverse waves, respectively. The curves are averaged 
over different directions of q, taking advantage of the 
the isotropy of the system. The error bars contain 
contributions from the finite optical resolution of the 
microscopy technique (assumed to be « 100 nm), which 
affects the amplitudes of modes in Fourier space. As 
expected, they show an initial quadratic growth (in 
agreement with Debye behavior) and exhibit a maximum 
near q = tt where the system is stiffest and propagation 
of waves fades out. 

Usually, the growing amplitude of the curves with 
increasing T reflects the stiffening of the system upon 
cooling. In a harmonic crystal, the “dispersion relation” 
scales directly with the coupling parameter since the 
spring constants grow quadratically with and thus go 
linear with T |22|. Accordingly, the curves for different 
temperatures should collapse onto one single master 
curve if divided by the interaction strength. Interest¬ 
ingly, we do not find this scaling in the amorphous 
solid. Besides the magnetic field, T only depends on 
intrinsically fixed parameters as laboratory temperature, 
magnetic susceptibility or area density. Since the 
particle-particle interactions (and therefore the “spring 
constants”) depend crucially on the respective distance, 
a subtle change in structure might explain this behavior. 
With increasing interaction parameter, the particles try 
to maximize their distances at constant Voronoi area in 
order to minimize energy. 

Because in the disordered state the equilibrium particle 
positions are not fixed for all times, the structure of the 
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FIG. 2. Dispersion relations of longitudinal and transverse ex¬ 
citations rescaled with the linear dependence on F expected in 
harmonic systems (hlled and empty symbols in top and bot¬ 
tom plot, respectively). The curves for different temperatures 
do not collapse, which can only be explained if structural re¬ 
arrangement occurs. Due to the better optical resolution only 
trajectories of big particles are analyzed here and in the fol¬ 
lowing plots. 

system might adapt to the increased pressure without 
affecting the macroscopic area density. This effect is 
barely visible in the structure factor. 

On a side note, the growth of the dispersion amplitude 
mirroring the increasing rigidity tackles another unique 
feature of our model system. Since the pair potential 
is purely repulsive and the system is investigated at 
constant volume (and particle number), the ensemble is 
intrinsically under pressure within the boundaries given 
by the sample cell. Since pressure equals an energy 
density (that is, energy per area in 2D), the pressure is 
given by the control parameter T, too, such that tem¬ 
perature and pressure can not be varied independently. 
A quantitative analysis (see Fig. shows that the value 
of the dimensionless pressure given by F is of the same 
order of magnitude compared to dimensionless elastic 
moduli. Such pressure affects the sound velocities given 
by the slope of the dispersion relation. 

Now, we turn to the main results, namely the elastic 
constants derived from the “dispersion relation”. Fig. 
shows exemplary the dispersion relation for F = 64 and 
F = 287 divided by according Eqn. and The 
solid curves are fits to the data to get the elasticity 
indicated by green arrows. We chose an intermediate 
regime (0.8 < qa < 2.5) for the extrapolation (f —>■ 0 



FIG. 3. jj, and 2^ + X are given by the intercept of an 
extrapolation (indicated by green arrows) derived from the 
longitudinal and transverse dispersion bands (filled and empty 
black circles) for a solid (blue) and fluid (red) sample. The 
gray shaded region is used for fitting the data and the inset 
shows the same plot with logarithmic scaling. 

in the shaded region where the data fits best. While 
in a linear fit was used to get the intercept we here 
expand ^ = 1 — ^ + 0{q‘^) and fit quadratically. 

30 measurement runs were performed at various tem¬ 
peratures. Due to prior work the transition was 
expected at 180 < T < 240. Therefore, we focused on 
good temperature resolution in that temperature range, 
resulting in high quality data each with measurement 
length of at least 8 • 10"^ seconds. The result is shown in 
Fig. i The error bars contain contributions from the 
finite optical resolution and the extrapolation of reduced 
dispersion relations towards small q as described by 
Eqn. and [^. The statistics for the latter was greatly 
improved by averaging over different directions of q. 
In the limits of these errors, the shear modulus below 
F « 190 equals zero, regardless of the applied external 
magnetic field. For increasing temperatures above, the 
moduli grow linearly, reflecting the linear “stiffening” 
of the interaction potential. With the amount of data, 
a zoom in on the shear modulus (shown in the lower 
plot) clearly reveals an abrupt rise of /i, consolidating 
the results of simulation m and delivering the first 
experimental evidence for a discontinuous behavior 
of the shear modulus at the glass transition. This 
jump separates the fluid phase from the amorphous 
state, as indicated by the shading from red to blue. 
Averaging about the data in the transition region, 
we find Ft = 195 ± 5. The scattering of the data is 
attributed to dynamical heterogeneity present in our 
system 1371 which have been shown to be correlated 
with elastic heterogeneity [38] . The dashed lines are fits 
to the moduli in the region 220 < F < 450. For the 
amorphous system we find a slope of /r ~ 0.31 • F and 
B = A -I- /i ~ 1.21 • F. This compares to measurements 
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FIG. 4. Temperature dependent bulk- and shear modulus. 
The data suggest a jump of moduli at the transition temper¬ 
ature, which is located at Ft = 195±5. The lower plot zooms 
in on the ordinate. In the solid, the moduli grow linearly due 
to increased pressure indicated by dashed lines and we extrap¬ 
olate this linear behaviour down to the fluid. The inset shows 
a zoom of ordinate and abscissa: the data are not compatible 
with a continuous increase of elasticity as e.g. given by the 
dotted lines. 


for dipolar monodisperse crystals^ where we found 
/i = 0.35 • r and B = 3.46 • F with a melting temperature 
at r = 60, see Ref. [S]. The latter is consistent with 
a thermodynamic calculation for a dipolar crystal at 
T = 0 which yields a ratio of bulk to shear modulus of 
B/^i = 10 [35] • An increased ratio B/yi compared to 
crystals is proposed by nonafline displacements affecting 
mainly the shear modulus of systems with short range 
interaction EOj. For the dipolar monolayer the ratio 
B/fi 3.5 is decreased, dominantly due to a reduced 
bulk modulus. 

Concerning the height of the jumps we are not aware of 
predictions of universal behavior of elasticity in thermal 
amorphous systems as e.g. the famous 167r « 50.26 for 
Youngs modulus in 2D melting [41] . but the measured 
values in glass and crystal are of the same order. For the 
jump in the bulk modulus BgoUd = w • Bfi^id a factor 
m = 1.7 is calculated by mode coupling theory for hard 
spheres in 3D [12] and implicitly given in [33] where it 
is shown that 3D and 2D systems behave comparably. 
For the dipolar system at hand we observe a factor 
TO « 6 for the glass and to « 4 for a crystal [H]. In 
comparison to hard discs we attribute this increased 
factor to the softness of the dipolar potential. This 



FIG. 5. Finite time effects investigated by varying the time 
interval At used to determine the particle “equilibrium posi¬ 
tion”. On the glassy side, the shear modulus does not decay 
to zero even for the longest integration times. Data of the 
system at F = 287 in the solid (blue) and F = 103 in the fluid 
(red) while the bulk and shear moduli are marked by filled 
and empty symbols. The inset shows the high frequency re¬ 
gion in the frozen state, which suggests a power law behavior. 
Upper and lower black lines are functions with exponents — | 
and — 


discontinuous behaviour distinguishes the glass transi¬ 
tion from jamming phenomena, for which the modulus 
was found to grow continuously from zero. This is 
encompassing granular systems [Hl^, percolating gels 
with attractive interactions [531 154] and furthermore 
systems exhibiting friction between particles m ES]. 
This might be an ubiquitous difference between thermal 
and granular systems. 

While a recent discussion about 2D and 3D glass claims 
to find fundamental differences [45] we expect the qual¬ 
itative behavior to be robust against changes in dimen¬ 
sionality. In 2D crystals the existence of Peierls insta¬ 
bilities [46] or Mermin-Wagner fluctuations [471 SH] as 
long wavelength density fluctuations (of the order of sys¬ 
tem size) is well established. Peierls argument using rel¬ 
ative distance fluctuations is not based on periodicity 
(just a typical particle distance is needed unlike e.g. in 
a gas) and thus we expect such long wavelength density 
fluctuations (with diverging amplitude in the thermody¬ 
namic limit) also in 2D glass. Using local coordinates 
namely cage-relative displacements (called modified dy¬ 
namic lindemann-parameter for 2D crystals [311 EH]) as 
suggested in [miE] would solve most of the discrepancies 
found in [45]. In the present analysis global coordinates 
are used and Mermin-Wagner fluctuations are expected 
to suppress the elastic response in the limit of infinite 
time scales. Thus the time interval At to determine the 
quasi-equilibrium position has to be chosen carefully. 

In Fig. [^ the finite time effects of the evaluation 
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FIG. 6. Graph a shows the spectrum of eigenvalues of the displacement field covariance matrix of big particles. By scaling 
with LJ, the inset reveals convincing Debye behavior. The red arrows indexed c-f indicate positions in the spectrum where the 
mode structure is shown in graphs c-f. The corresponding “frequencies” are given above the respective graph. For details, see 
text. 


procedure are shown. By changing the averaging time 
At determining the quasi-equilibrium positions, different 
frequency regimes are probed. Short averaging times 
yield a significant increase in elastic moduli, reflecting 
the expected stiffening at high frequencies even for fluids 
[3]. For low frequencies (large At), in contrast, the shear 
modulus of the fluid state decays to zero, while in the 
glassy state, it decays into a plateau, even for times 
up to 10® seconds. This stability of the shear modulus 
is remarkable, because other quantities measuring the 
system dynamics already show signs of relaxation. This 
averaging time used in Fig. |4 is At = 1.8 • 10“^ indicated 
by the green arrow in Fig. ^ (limited by the longest 
accessible time from the fluid data). Comparing with 
Fig.g this is well after the onset of the alpha relaxation 
time Ta for the datasets in the vicinity of the transition. 
In contrast to the shear modulus, the bulk modulus is 
much more sensitive to increased dynamics, as it does 
not show a well developed plateau for low frequencies. 
This is attributed to the fact that transverse shear waves 
have much larger amplitudes than longitudinal compres¬ 
sion waves, such that distortions of the wavefront affect 
properties derived from the latter much more strongly. 
This striking behavior proves the shear modulus to be 
a good quantity to investigate the glass transition. The 
inset finally shows a log-log plot of the high frequency 
regime for the frozen state. For about two decades, 
both moduli show a power law behavior. The blue lines 
are functions with exponents of — describing high 


frequency bulk and shear properties, respectively. This 
power law behavior challenges theoretical descriptions. 

In the last section we turn to the analysis of the shape 
of the excitations. Fig. |6]a) shows the spectrum of eigen¬ 
values D(u}) of the displacement field covariance matrix 
Cij for the system at F = 200, a temperature slightly be¬ 
low the glass transition temperature. The most principal 
components (carrying the most weight) correspond to low 
“frequencies” due to the mapping uj oc Because 

Cij has the dimensions of square meters as it measures 
amplitudes, the “frequencies” have the unit Only 

a valid mapping of amplitudes, wavelength, and frequen¬ 
cies (implicitly using a dispersion relation) connects a 
length scale given by the eigenvalues of the covariance 
matrix to frequencies used in the density of states D{uj). 
An interesting feature of D{oj) is the double peak struc¬ 
ture, centering at a; « 0.8 p.m~^ and w « 2.1 pm~^. 
We interpret this as a reminiscence of van-Hoove type 
singularities as they appear e.g. in crystals due to the 
vanishing sound velocity at the Brillouin zone edge [35]. 
Naively one could argue that first and second peak repre¬ 
sent longitudinal and transverse modes respectively, but 
since PCA is a priori insensitive to polarizations the con¬ 
nection is less trivial. Nevertheless, Fig. I^b) shows the 
density of states as function of temperature. The peaks 
shift to lower frequencies with increasing temperature 
(decreasing interaction strength) reflecting a softening of 
the system as expected. Note that the first peak grows 





















7 


at the expense of the second upon approaching the fluid 
phase. This resemblance with a density of states known 
from a crystalline solid also shows up in the reduced den¬ 
sity of states D{uj)/uj which is shown in the inset of Fig. 
I^a). Scaling with w yields a constant value for the low 
frequency region up to the first peak. This is in accor¬ 
dance with the Debye model D{uj) oc , where d is the 
dimensionality of the investigated system. In the Debye 
model, it is plane wave excitations which make up the low 
frequency part of the spectrum. Putting issues regarding 
the amplitudes frequency mapping aside, the mode struc¬ 
ture is clearly resolved by the eigenvectors of the most 
principal components of the displacement field. For de¬ 
creasing wavelength, these are shown in graphs |^c-f with 
the corresponding frequency noted above. Those frequen¬ 
cies are also indicated by red arrows in Fig. §a). While 
|^e,d) show some kind of Lissajous-pattern, for modes 
of a given frequency but different phases, they start to 
show incoherence for higher frequencies indicated by a 
random distribution of arrows. The coherence length of 
the imaged virtual motion gets notably smaller, reflect¬ 
ing the growing influence of disorder. This culminates 
at frequencies accounting to the first peak, where the 
mode structure appears almost disordered in Fig. §f). 
Nonetheless, care has to be taken since PCA was shown 
to generate mode pattern from mixed states at higher 
frequencies; the low frequency modes on the other hand 
are expected to be pure states [SS]. In Fig. |^c) one 
can clearly verify the plane wave character of the excited 
states in the low frequency limit. Accordingly, the ansatz 
to equally distribute energies via the equipartition theo¬ 
rem on plane wave like modes for lim ((7 —>■ 0) (Eqn. is 
reasonable. 


CONCLUSION 

We have experimentally investigated the elastic prop¬ 
erties of a colloidal glass former in great detail. By tun¬ 
ing the interactions in situ, we have accessed the well 
defined displacement field in Fourier space at numerous 
effective temperatures and used the equipartition theo¬ 
rem to calculate the dispersion relations. Scaling with the 
interaction parameter reveals a hidden structural change 
upon cooling which is too subtle to be picked up by other 
functions. In the limit of long wavelength, the disper¬ 
sion relations yield the elastic moduli of compression and 
shear, of which particularly the latter is of great interest 
for an deeper understanding of the glass transition. It 
has been demonstrated that the shear response unam¬ 
biguously separates the fluid state from the amorphous 
solid. The transition itself exhibits a discontinuous be¬ 
havior, which is in agreement with theoretical predic¬ 
tions and simulations and distinguishes it from granular 
matter or gelation type scenarios. Finite time analysis 
reveals a striking stability of the shear-related measure¬ 


ments, which are proven to be decoupled from relaxation 
phenomena like cage breaking. This makes the shear 
modulus an exceptionally good quantity to investigate 
the glass transition. 
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